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Abstract 



It has been known for some time that proportional output feedback 
will stabilize MIMO, minimum-phase, linear time-invariant systems if the 
feedback gain is sufficiently large. High-gain adaptive controllers achieve 
stability by automatically driving up the feedback gain monotonically. 
More recently, it was demonstrated that sample-and-hold implementations 
of the high-gain adaptive controller also require adaptation of the sampling 
rate. In this paper, we use recent advances in the mathematical field of 
dynamic equations on time scales to unify and generalize the discrete and 
continuous versions of the high-gain adaptive controller. We prove the 
stability of high-gain adaptive controllers on a wide class of time scales. 
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1 Introduction 

The concept of high-gain adaptive feedback arose from a desire to stabilize cer- 
tain classes of linear continuous systems without the need to explicitly identify 
the unknown system parameters. This type of adaptive controller does not iden- 
tify system parameters at all, but rather adapts the feedback gain itself in order 
to regulate the system. A number of papers examine the details of various kinds 
of high-gain adaptive controllers [TTJ [15l HH [23] , among others. More recently, 
several papers have discussed one particularly practical angle on the high-gain 
adaptive controller, namely how to cope with input/output sampling. In par- 
ticular, Owens [17] showed that it is not generally possible to stabilize a linear 
system with adaptive high-gain feedback under uniform sampling. Thus, Owens, 
et. al., develop a mechanism to adapt the sampling rate as well as the gain, a 
notion subsequently improved upon by Ilchmann and Townley [121 1131 114j , and 
Logemann |16j . 
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In this paper, we employ results from the burgeoning new field of mathe- 
matics called dynamic equations on time scales to accomplish three principal 
objectives. First, we use time scales to unify the continuous and discrete ver- 
sions of the high-gain controller, which have previously been treated separately. 
Next we give an upper bound on the system graininess to guarantee stabilizabil- 
ity for a much wider class of time scales than previously known, including mixed 
continuous/discrete time scales. Third, the paper represents the first applica- 
tion of several very recent advances in stability theory and Lyapunov theory for 
systems on time scales, and two new lemmas are presented in that vein. We 
also give a simulation of a high-gain controller on a mixed time scale. 



2 Background 

We first state two assumptions that are required in the subsequent text. 

(Al) The system model and feedback law are given by the linear, time- invariant, 
minimum phase system 

x(t) = Ax(t) + Bu(t), y(t)=Cx(t), x(0) = x , (1) 
u(t) = -k(t)y(t), (2) 

for t > 0. System parameters A £ R" xn , B,C T G R nxm , x £ R" and 
n are unknown. The feedback gain k : R — > R + is piecewise continuous, 
and nondecreasing as t — > oo. By minimum phase, we mean that the 
polynomial 

\A - XI B 
C 



det 



with A £ C is Hurwitz (zeros in open left hand plane). 
(A2) Furthermore, 

{CB) T + (CB) > 0, (3) 

i.e. (CB) T + (CB) is positive definite. (In [TB] it is pointed out that a 
nonsingular input/output transformation T always exists such that B = 
BT- 1 and C = TC give [CB) T + (CB) > 0.) 

Under these conditions it has been known for some time (e.g. [TT]) that 
there are a wide class of gain adaptation laws k{t) — f(y(t)), f : R m — * R, that 
can asymptotically stabilize system ((TJ) in the sense that 



y £ L 2 [io,oo), lim k(t) < 



t — >oo 



■yj. 



Subsequently, various authors [TTJ [12l [T71 [19] assumed that the output is ob- 
tained via sample-and-hold, i.e. yi := y(ti) with k% := fe(ij) and i £ No- Thus 
it becomes necessary also to adapt the sample period hi := U + \ — U so the 
closed-loop control objectives are 

yi £ £ 2 {io,oo), lim fcj < oo, lim hi > 0. 
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Though several variations on these results exist, these remain the basic control 
results for continuous and discrete high-gain adaptive controllers. The contin- 
uous and discrete cases have previously been treated quite differently, but we 
now construct a common framework for both using time scale theory. 



3 A Time Scale Model 

The system of (Al) can be replaced by 

x A (t) =A(t)x(t) + B{t)u(t), y{t)=Cx(t), x(0)=x , (4) 
u(t) = k(t)y{t), t£T, (5) 

where T is any time scale unbounded above with € T. With a series expansion 
similar to pj] , we see that 

A(t) := expc(fi(t)A)A, B(t) = expert) A) B, (6) 

where expc is the matrix power series function 

expc(X) := / + -X + \ ; X 2 + • • • + -X^ 1 + ■■■ (7) 
2 6 n\ 

and fi is the time scale graininess. Implementing control law ([2]) then gives 

x A (t) = A(t)x(t), A{t) := expc{n(t)A)(A - k{t)BC). (8) 

Note that A, B and A may all be time- varying, but we will henceforth drop 
the explicit reference to t for these variables. For future reference, we also note 
that if sup tgT /i(t) < oo, then ||j4|| and ||B|| are also bounded (c.f. Appendix, 
Lemma 7.4). 

The design objectives are to find graininess ix and feedback gain k as func- 
tions of the output y, 

H(t) = g(y(t)), 5 :M m ^[0,co), (9) 
fc(t) = f(y(t)), f : R m -» M+, (10) 

with < 6 T and k £ C r( j(T) nondecreasing, such that 

y £ L 2 [t , oo) T := |w : \w(t)\ 2 At < ooj , Urn k(t) < oo. 

It is important to keep in mind the generality of the expressions above. A 
great deal of mathematical machinery supports the existence of delta derivatives 
on arbitrary times scales, as well as the existence and characteristics of solutions 
to ([5]). See, for example, e.g. 
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4 Stability Preliminaries 

We begin this section with a definition and theorem from the work of Potzsche, 
Siegmund, and Wirth [20] : 

Definition 4.1. The set of exponential stability for the time-varying scalar 
equation z A (t) = X(t)z{t), zq :— z(to) € C with A : T — > C and z,Ae C r d(T), is 
given by 

5(T) := 5fc(T) U 5 E (T) 



re 



S C (T) := U) G C : a = - limsup -J— /* i^^±^MM A< > 

(11) 

5 R (T) := {r/(i) G R : Vi G T, 3r >< uni/i r G T suc/i i/iat 1 + m(*M*) = 0} , 
loii/i n G C r( j(T) and arbitrary ijET. 

Theorem 4.2. |20j Solutions of the scalar equation z A (t) = A(i)z(t) are expo- 
nentially stable on an arbitrary T z/ and onZj/ if X(t) G <S(T). 

We note here that Potzche, Siegmund, and Wirth did not explicitly consider 
scenarios where r\ is time- varying, but their stability analysis remains unchanged 
for rj{t). 

The set 5r(T) contains nonregressive eigenvalues A(i), and a loose interpre- 
tation of 5c (T) suggests that it is necessary for a regressive eigenvalue to reside 
in the area of the complex plane where ||1 + (i(t)X(t)\\ < 1 "most" of the time. 
The contour ||1 + fj,(t)X(t)\ | = 1 is termed the Hilger Circle. Since the solution 
of z A {t) = X(t)z(t) is x(t) = z e\(t, t ), Theorem 4.2 states that, if A(t) G 5(T) 
then some K(to) > 1 exists such that 

\\z(t)\\ = \\z \\\\e x (t,to)\\<\\z \\Ke- a ^ t °\ (12) 

where e\ is a generalized time scale exponential. The Hilger Circle will be 
important in the upcoming Lyapunov analysis, as will the following lemmas. 

Lemma 4.3. Let z G C r d(T) be a function which is known to satisfy the in- 
equality z A (t) < \{t)z(t), 2(0) = z Q G K with A : T -> K and A G C rd (T). 7/ 
z(<) > for all t>t eT, then X(t) G 

Proof. Defining u(f) = z A (t) — X(t)z(t) gives rise to the initial value problem 
z A {t) - X(t)z(t) = v(t), z(Q)=za, teT, (13) 

where v(t) < 0. 

First, suppose A(t) is negatively regressive for t > t , i.e. X(t) < —l//i(t) 
with fj,(t) > 0. Then ||T3J| yields z a = fiXz + z + fxv = z(l + fj,X) + \iv < 0. 
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On the other hand, suppose A(T) is nonregressive for some T > to, T G T. 
If A(i) < — l/n(t) over T > t > to, then invoke the preceding argument. If 
A(t) > over T > t>t , then solve (H3J) to get 

z(t) = e x (t,to)z + [ e x {t,a(T))v(T)AT, T>t>t . (14) 

Since eA(i,io) > for t < T, we see that J t eA(i, g(t))v(t)At < 0. However, 
for t > T, HH becomes 

z(t) = + e x (t,a(T))v(T)AT, t>T (15) 

J to 

Thus, z(t) < for all t > t for both negatively regressive and nonregressive X(t), 
a contradiction of the Lemma's presupposition. This leaves only X(t) G 1Z + . □ 

At this point we pause briefly to discuss Lyapunov theory on time scales. 
DaCunha produced two pivotal works [5J [5] on solutions P of the generalized 
time scale Lyapunov equation, 

A(t) T P(t) + P{t)A{t) + fi(t)A T (t)P(t)A(t) = -Q(t), t G T, (16) 

where A(t), P(t),Q(t) G R nxn , A and Q are known and Q(t) > 0. Though it 
will not be necessary to solve (|16|) in this worlQ, we will see that the form of 
(|16p leads to an upper bound on the graininess that is generally applicable to 
MIMO systems, an advancement beyond previous works which gave an explicit 
bound only for SISO systems. 

Before the next lemma, we define 

JC := {k : T -> R+, ke C rd {T), k A (t) > Vi G T, lim k(t) = oo}. 

t — >oo 

The next lemma follows directly. 

Lemma 4.4. Given assumptions (Al) and (A2) and fc G /C, i/tere exists a 
nonzero graininess £t(t) and a time t* such that, for all fi(t) < /2(£) and t > t* , 
the matrix —kCB satisfies a time scale Lyapunov equation with P — L , Q(t) > 
£2 for small 62 > 0, and B from @ . 

1 DaCunha proves that a positive definite solution P(t) to the time scale Lyapunov equation 
with positive definite Q(t) exists if and only if the eigenvalues of A(t) are in the Hilger circle 
for all ( 6 T. Furthermore, P(t) is unique. As with the well known result from continuous 
system theory (c.f. |210 . the solution is constructive, with 

P(t) = f >S> AT(t) (s,0)Q(t)<£ A{t) (s,0)As, 
J s £ 

where <E>A(i, to) denotes the transition matrix for the linear system x^(t) = A(t)x(t), x(t$) = 
I. The correct interpretation of this integral is crucial: for each t E T, the time scale over 
which the integration is performed is St := fi(t)No, which has constant graininess for each 
fixed t. 
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Proof. We construct /i(t) as 

_ m = _±_ ( \ min {CB + (CB) T } 
M ' ' k(t) \ \ max {(CB) T CB} 

with E\ > sufficiently small so that > on T. This holds if 

(CB) T + CB- fl(t)k{t)(CB) T CB = sil. (18) 

Multiplying HU) by -k(t) gives 

(-k(t)CB) T + (-k(t)CB) + fd(t)(-k(t)CB) T \-k(t)CB) = -k(t)eil. 

(We now drop the explicit time-dependence for readability.) Set 

E(/2) := expc(flA) — I, 

so that B(p) = [I+ S(/2)]B, yielding 

(-kCB) T + (-kCB) + fi(-kCB) T (-kCB) = -k{ej - Z), 

where each term of Z is a product of contants times S(/2). Since | |E(/2)| | — > as 
fi — > (c.f. Appendix, Lemma 7.4), there exists a time t* 6 T when \\Z\\ < E\. 
Because the preceding arguments admit any graininess < p(t), it follows 
that 

(-kCB) T + (-kCB) + n(-kCB) T (-kCB) = -kQ, t > t*, fj,(t) < p,(t), 

(19) 

where Q(t) > e 2 and e 2 := £i - ||Z|| > 0. □ 

We comment briefly on the intuitive implication of Lemma 4.4. Equation 
(flT)|) shows that there exist positive definite Q(t) and P = I to satisfy an equa- 
tion of the form (p~6|) . According to DaCunha, this implies that the eigenvalues 
of —kCB lie strictly within the Hilger circle for t > t* . 



t e T, 



(17) 



5 System Stability 

We now come to the three central theorems of the paper. If BC is not known 
to be full rank, or cannot be full rank because of the input/output dimensions, 
then it must be assumed (or determined a priori) that the eigenvalues of A — 
k{t)BC attain negative real parts at some point in time. This phenomenon 
is investigated in-depth by other authors [TTJ [53]. We then make use of the 
observation by Owens, et. al. [18] . that there must exist some k* > such that 
if k[t) = k* the system of (Al) has a positive-real realization. This, together 
with the Kalman-Yakubovich Lemma [21] , implies existence of P, Q > so that 

(A - k*BC) T P + P(A — k*BC) — —Q and PB = C T CB. (20) 
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Theorem 5.1 (Exponential Stability). In addition to (Al) and (A2), suppose 

(i) t £ T where T is a time scale which is unbounded above but with /i(t) < 

(ii) k £ K, [implying from (| 1 T[) that fj,(t) — ► 0, 6u£ not necessarily monotoni- 
cally), 

(iii) PC is not necessarily full rank, but there exists a time t* £ T such that 
the eigenvalues of (A — k(t)BC) are strictly in the left-hand complex plane 
for t>t*. 

Then the system ijl j) . ([5 ]) is exponentially stable in the sense that there exists 
time to £ T and constants K(to) > 1, a > 0, such that 

\\x(t)\\ < 11^0)11^-^*-*°), t>t Q . 

Proof. Set k* := k(t*). Then assumption (iii) is the prerequisite for equation 
(|20|) . Again, we suppress the time-dependence of x, k, /j,, A and B. Similarly 
to Lemma 4.4, terms containing £(//) may be added to the first equality in (|20p 
to obtain 

(A- k*BC) T P + P(A- k*BC) < -£ 3 t>t*\ (21) 

for some small £3 > 0. Note t — t** is the point at which terms involving 
become small enough for (JUJ) to hold. Defining Z x := C T C'Y.B - PEP, the 
second equality in (|2H)) gives 

PB = CFCB + Z-l. 

Consider the Lyapunov function V — x(t) T Px(t) with P from (f2T)]) . Then, 
using [pO]). HH and Lemma 4.4, 

= x T [yl T p + p»4 + M T -P^ 

= x T \(A - k*BC) T P + P(A - k*BC) -(k- k*){C T B T P + PBC) 
+H(A - kBC) T P{A - kBC)] x 

< x T [-e 3 - (k - k*)(C T B T C T C + C T CBC + 2Z X C) 

+ [iA T PA - 2[ikA T PBC + nk 2 C T B T C T CBC + ^k 2 C T B J ' Z x C\x 
= x T [~e 3 + fiA T PA}x 

+ y T [-k{CB) T - kCB + fik 2 (CB) T CB + fik 2 B T Z 1 + k*((CB) T + CB)]y 
+ x T [-(k - k*)2Z x - 2[ikA T PB\y 

< -(e 3 -fi\\A T PA\\)x T x 

-(k£ 2 -fik\\B\\kZ 1 \\-2k*\\CB\\)y T y 
+ (2\\kZ 1 \\ + 2^k\\A T PB\\)x T y. 
At this point we observe the following: 
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• ||A T PA|| is bounded because A is bounded. Set 71 := sup teT ||A T PA||. 

• fik is bounded by assumption (i) and (jTTJ) . 

• ||fcZi|| is bounded because k is proportional to — and — ► const 
as /1 — > 0. Set 72 := sup tgT 2||feZi||. 

• Set 73 := sup teT /ifc||S||7 2 . 

• Set 74 := sup teT 2fc*||C*i?||. 

• Set 75 := sup feT 2^fc||i T Pi?||. 

Recalling the standard inequality x T y < f3x T x + jjV T y for any (3 > 0, we 
then have 

V A < -(£3 - M7i - Pl2 - (3l5)x T x - (ke 2 - 73 - 74 - ^72 - ^lh)y T y- 

By assumption (ii), there exists a time to > max{i*,£**} such that, for suffi- 
ciently small (3, 

yA < -( g 3-M7i- g7 2 -/?75V := vm t > ^ 

with r/(t) < 0. Then, by [2] Theorem 6.1], Theorem 4.2, and Lemma 4.3 it 
follows that there exists K(to) > 1 so that 

||*(i)|| <if||^o)l|e-^ (t - to) , ^><o, 

where a is defined in (jTTJ) . □ 

We point out that, when BC is full rank, assumption (iii) above is no longer 
necessary as there always exists a k* such that the eigenvalues of (A — k(t)BC) 
are strictly real-negative for k(t) > k* . One more lemma is required before the 
next theorem. 

Lemma 5.2. IfT is a time scale with bounded graininess (i.e. /ioo := sup tgT /j,(t) < 
00 ), then 

ci / e at dt< e at At<c 2 e at dt, 

J to J tg J tQ 

where c\, c 2 , a £ R and C\, C2 > 0. 

Proof. Consider the case when a > 0. The process of time scale integration 
is akin to the approximation of a continuous integral via a left-endpoint sum 
of (variable width) rectangles. If the function to be summed is increasing (as 
in this case), the sum of rectangular areas will be less than the continuous 
integral, meaning c 2 — 1, c\ < 1. One estimate of the lower bound, then, follows 
by simply increasing c± until c\e at meets one of the rectangle right endpoints 
which are given by e ap ^\ Thus c\e at < e ap ^\ or equivalently, ae " 7 ^' < e at . 
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This in turn yields c\ < e - "^*'. Therefore, the most conservative bound is 
given by c\ < e~ aMo °. The case for a < can be argued similarly, leading to 
the lemma's conclusion: 

(e- af *°°, a>0, fl, a>0, 

ci = < c 2 = < 

[1, a<0, [e~ Q ^, a < 0. 

□ 

We are now in a position to state the main theorem of the paper. 

Theorem 5.3. In addition to (Al) and (A2), assume the prototypical update 
law, k A (t) — \\y(t)\\ 2 with ko := fc(0) > 0. Then lim^oo k(t) < oo and y € 
L 2 [t , oo)t- 

Proof. For the sake of contradiction, assume k(t) — > oo as t — > oo. Then k E K,. 
Theorem 5.1 yields x € £oo[^o> oo)t and therefore y € ioo[^o, oo)t. The solution 
for k(t) is (by Theorem 2.77]), 



fc(t) - k(t ) + f \\y(T)\\ 2 Ar < k(t ) + \\C\\ 2 f ||x(r)|| 2 Ar, 

Jtn Jtn 



t > t 



0. 



'to 

In conjunction with Lemma 5.2, this allows 



lim kit) < k(t ) + \\C\\ 2 / \\x{t)\\ 2 At 

f - > °° Jto 

<k(t ) + \\C\\ 2 \\x(t )\\ 2 K 2 / e-^-^At 



= k Q + \\C\\ 2 \\x(t Q )\\ 2 K 2 
< oo. 

This contradicts the assumption, so it must be that k(t) < oo for t 6 T. It also 
immediately follows that J | \y\ \ 2 At < oo. □ 

It seems possible that Theorem 5.3 may be improved to show that the output 
is convergent, i.e. y(t) — > as t — > oo. This is left as an open problem. 



6 Discussion 

We remark here that there is a great amount of freedom in the choice of the 
update law for k(t) (c.f. [E]). We use the simplest choice for convenience (as 
do most authors); the essential arguments remain unchanged for other choices. 
There is also freedom in the choice of update for fi(t). The expression (fT5)l . can 
be simplified to 

(CB) T + CB-n{t)k{t)(CB) T CB>0, teT, 
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for any fi(t) < p,(t). In the SISO case, this further reduces to the expres- 
sion derived by Owens [17], that [ikCB < 2. It requires the graininess (which 
may interpreted as the system sampling step size for sample-and-hold systems) 
to share at least an inverse relationship with the gain, but is otherwise quite 
unrestrictive. Ilchmann and Townley [T3] note that fi(t) — m \^%m meets 
fikCB < 2 after sufficient time without knowledge of CB. While previous 
works have always constructed a monotonically decreasing step size, the time 
scale-based arguments in this paper reveal even greater freedom: /i(t) may ac- 
tually increase, jump between continuous and discrete intervals, or even exhibit 
bounded randomness. Two examples of the usefulness of this freedom come 
next. 

For the first example, we point out that the notation in the previous sections 
somewhat belies the fact that the system's time domain (its time scale) may 
be fully or partially discrete, and thus there is no guarantee in Theorem 5.3 
that the output has stabilized between samples. As pointed out in [13] and 
elsewhere, a sampled system with period [i is detectable if and only if 

^r^/i i Z for any A fc ^ A*; j = V=I, (22) 
2ttj 

where Afc,A; £ {0,spec(A)}. We next comment on how to circumvent the in- 
trasample stabilization problem. Both of the following methods essentially per- 
mit the graininess fi to "wiggle" a bit so that an output sample must eventually 
occur away from a zero-crossing. Recalling that fi(t) < p,(t) from Theorem 
5.1 (i), let [i(t) = p,(t)v(t), where < v < 1 is one of the sequences below: 

1. ) Let v be consist of an infinitely repeated subsequence with n! + 1 elements 
that are random numbers between and 1. Let these elements, labeled Vi, V2, 
be irrational multiples of each other. Assume [i(t) converges to fivi > such that 
y(t) = but the true continuous output is nonzero. This implies that there exist 
integers k,l such that (Afc ~^ )p '" 1 S Z. As sequence v advances, there may be at 
worst nl combinations of j, k such that *- Aj ^j p ' Vr <E 7L for r = l...n\. However, 
at the next instant in time, ( Afc ~^ ) '^ t '" l + 1 must be irrational and therefore not 

' 27TJ 

in Z. The controller will detect a nonzero output and continue to adapt k(t) 
and /i(t). In practice, of course, it is not possible to obtain a sequence of 
truly irrational numbers, but most modern computer controllers have enough 
accuracy to represent the ratio of two very large integers, so that this technique 
would only fail for unpractically high magnitudes of ||Afe — Aj||. 

2. ) Let v be a sequence of random numbers in the specified range. Even in 
a computer with only 8-bit resolution for v, the probability of € Z 
drops drastically after a few sample periods. 

We remark here that, if (|22|) holds then (A, C) is detectable because (A, C) 
is detectable. Thus, the stability of y implies the stability of x. We do not dwell 
on this here, but see e.g. [13] for a similar argument. 

For the second example, we consider a problem posed by distributed con- 
trol networks (c.f. [U [5]). Here, one communication network supports many 
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ft(t) 
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Figure 1: In graph (A), output y(t) (thick line) and gain k(t) are plotted. In 
graph (B), [i(t) (thick line) is kept at 90% of jl{t) during blocking. A block 
occurs after every 1 second of continuous control. 



control loops as well as a certain volume of unrelated high-priority traffic. The 
(unpredictable) high-priority traffic may block the control traffic at times, forc- 
ing longer-than-anticipated sample periods. In normal operation, the controller 
may sample fast enough to behave, for all practical purposes, like a continuous 
control with fi(t) = 0. At the blocking instant tb, fi(tb) rises to some unpre- 
dictable level. The scheduling question is, when should the blocked controller 
emit a communications packet of high enough priority to override the block? 
One answer is straightforward: just before n(tb) exceeds fi(tb), the known max- 
imum sample period that will guarantee plant stability. (In [5] we suggest that 
even longer delays are possible under certain conditions.) Intuitively, lower gains 
permit longer blocking delays. 

The example of the previous paragraph is closely modeled using the variable 
P[a, b] time scale, which is continuous for an interval a, then has a gap for an 
interval b, then repeats. Figure [1] shows the regulation of a system implementing 
an adaptive gain controller in a blocking situation: 

C=[1,0}. (23) 

The gain begins at kg = ^ and the sampling period at /i(0) = 0. The bounding 
function for the graininess is fi(t) — ttJy (so that flkCB < 2). 



1 
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7 Conclusions 

In summary, the paper illustrates a new unified continuous/discrete approach 
to the high-gain adaptive controller. Using recent developments in the new 
field of time scale theory, the unified results reveal that this type of feedback 
control works well on a much wider variety of time scales than explored in 
previous literature, including those that switch between continuous (or nearly 
continuous) and discrete domains or those without monotonically decreasing 
graininess. Furthermore, several results relating to Lyapunov analysis on time 
scales appear here for the first time, including Lemma 4.3 (and its use in the 
proof of Theorem 5.1) and 5.2. A simulation of an adaptive controller on a mixed 
continuous / discrete time scale is also given. It is our hope that time scale theory 
may find wider application in the broad fields of signals and systems as it seems 
that many of the tools needed in those fields are beginning to appear in their 
generalized forms. 

We thank our colleague, Robert J. Marks II, for his very helpful suggestions 
throughout this project. 

8 Appendix 

We comment on the properties of the "expc" function referenced in the main 
body of the paper. 

Lemma 8.1. The power series ([7]) has the following properties: 

1. expc(X)X = Xexpc(X). 

2. expc(X) = (e x - TjX^ 1 when X^ 1 exists. 

3. For real, scalar arguments x, expc(x) = e 1 ^ sinc(^), where sine denotes 
the sine cardinal function. ( This is the motivation for the expc notation.) 



Proof. Parts 1-3 follow immediately from the definition. To verify 4, note 



Part 5 follows from a similar argument. Note that, by property 5, the decom- 
position expc(fiX) = I + gives ||S(/x)|| — > as /i — > 0, with uniform 



4. HexpcpOH <exp(||X||). 



5. ||expc(X) -7|| < exp(||X||) - 1. 




convergence. 



□ 
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